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Abstract 

In the discussion about the quantumness of NMR computation a con- 
clusion is done that computational states are separable and therefore can 
not be entangled. This conclusion is based on the assumption that the ini- 
tial density matrix of an individual molecule coincides with whole sample 
molecules distribution over single molecule energy levels. This means that 
quantum stochasticity is replaced by classical stochasticity. In the present 
paper it is shown, that quantum NMR computation can create genuine 
entangled states if initial system states are thermodynamical equilibrium 
ones. A separability analysis problem can arise when one interprets the 
readout signal from whole sample. 

PACS No.: 03.65.-w, 03.67.-a, 76.60.-k. 

Introduction 



It is commonly believed, that one of the reasons of the high effectiveness of 
quantum computations in comparison with the classical ones lays in the usage 
of the so- called entangled states of quantum objects. For example, entangled 
state of two quantum objects 

|a| 2 + |6| 2 = l. W 
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is constructed from the separate objects wave functions \ipi) and ^2} in such a 
way, that it can not be represented as a product of wave functions of the objects 
- that is it is non separable. 

Recently the works |l]]-[|3| have been published, from which it follows, that 
the density matrices of NMR system lay in the proximity of the maximally mixed 
state and that is why they are separable. As a result the NMR computations 
should be less effective comparing to what was expected before. In the present 
paper it is shown, that in suitable NMR experiments there is an entanglement 
in full extent and that the separability problem can arise during interpretation 
of the whole sample readout signal. 

2 The arguments in favor of NMR system states 
separability 

When analyzing the separability problem and possible entangled states absence 
in the NMR quantum computing usually the following reasoning is given. 

In the NMR quantum systems a cluster of a small number N of interact- 
ing spins (for convenience let us call it "molecule" ) plays a role of quantum 
processor. This molecule does not interact with other molecules of sample. 
The observed signal is a sum of signals from M molecules- "computers" of a 
whole sample. In liquid state NMR the values AE of typical intervals between 
molecule energy levels satisfy the condition AE << kT. Then the energy levels 
population differences are of the order of Me (where e — AE/kT « 1) and 
represent a small part of the total number of molecules M. At such conditions 
all of spin energy K = 2 N levels of individual processor are occupied approxi- 
mately equally and spin system initial state is described by the mixed density 
matrix 

Pm = EtlPk\<P k ra)W, ( 2 ) 
where Pk is the independent of m possibility to find a molecule in a pure spin 
state |</>m)- A quantum system state is called a maximally mixed state if all 
probabilities equal Pk = 1/K. Thus in liquid state NMR experiments the 
small parameter e means that during computations the system states are in the 
proximity of the maximally mixed state. 

In the paper |jj a mathematical proof is given that there is a small proximity 
of the maximally mixed state where density matrix is separable, that is it can 
be represented as a sum of products of individual particle density matrices. In 
the paper the lower limit of this proximity is given, whereas in the paper Q 
its upper limit is given too. These mathematical proofs are irreproachable. It 
means that the NMR quantum computation effectiveness is much less than in 
the case when states are entangled. 
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3 The arguments in favor of the NMR system 
states entanglement 



We think that the above mentioned separability analysis does not take into full 
account the underlying physics of the NMR experiments. The matter of fact is 
that during computations each processor carries on calculations independently 
of the states of other molecules. And only at the readout step the individual 
molecules signals are summed up. To interpret the results of quantum mechani- 
cal measurements as usually one must take into account that one deals here with 
"double stochasticity" - the quantum and the classical ones 0. The quantum 
one should be used when one averages over quantum object states which are 
defined by quantum wave functions. Whereas the classical one should be used 
when one averages over initial states of molecules. 

To adequately describe the NMR experiments one should perform two steps: 

a) to track down the detailed quantum dynamics of individual "processor" , 
each time starting with a possible initial state (using, for example, the wave 
function language); these steps should be repeated for all possible initial states; 

b) to average results over all possible initial states of all M molecules- 
"processors" (using density matrix formalism or an appropriate classical dis- 
tribution - the Boltzmann one, the Maxwell one and so on). 

Following this plan one can note, that at each time the state of an individual 
molecule in general is given by wave function 

l<U = £f=ia(Cbd>, (3) 

where a(t)^ n - numerical coefficients and \xm) i s the fc-th eigenstate of the 
molecule spin energy operator. 

It is well established in the NMR experiments fact that in the thermody- 
namical equilibrium a molecule occupies one of its stationary states 

\4>m) eq = \xt)i (4) 

which in the density matrix formalism can be written as p m ,eq = \Xm) (Xm\- 
This fact follows also from contemporary decoherence theory || , which describes 
how an open quantum system interacts with environment. According to the 
theory the coherent state (||) relaxes to one of the stationary state (^). 

During quantum computation such a molecule state may be transferred to 
non separable, entangled state by a few ways and so this state will be involved 
in computation independently of the states of other "processors" . These calcu- 
lations are the subject of quantum dynamics and not of statistical mechanics. 
The necessity in the latter arise when one reads out a result, which is a sum of 
individual molecule signals. 
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4 Signal from whole sample 



A sample state is defined by the direct product of individual molecule wave 
functions 

l*>eg = l<Meg ® l^)e, • ■ ■ ® = ®%=l \4>m) eq (5) 

Suppose that a sample is subjected to a sequence of external influences, 
which realize a quantum algorithm. The corresponding propagator is equal to 
Ucomp = Ilit/j, where C/j is a unitary transformation, which realizes a logic gate. 
And if one of these propagators transforms a molecule state to an entangled 
one, then the other gates will continue to work with such a state. 

Under the influence of propagator U comp the initial spin system state of the 
whole sample (^) transforms to 

I*}* = l<M' ® \H? ■ •• ® I^m)* = ®i f =1 |<U* (6) 

where = f7 COTOp |x„)- 

Suppose that to readout a result one has to measure the expectation value 
of some operator J x . Operators of physical quantities are additive 

^ = £™=i^, ( 7 ) 
and that is why in the basis (^) the expectation value of operator J x is 

_spM I k m \TJ+ Jxtj I k m \ \°) 

Since the sample molecules are supposed to be identical, then the right hand 
side matrix elements of Eq. (|^) depend only on the initial state index k m and 
not on the molecule number m. For all molecules with the same initial state k 
the matrix elements are the same and as a consequence the sum over m in the 
right side of Eq. (@) is reduced to the sum over k 

< J X >= 52fe=l Cfc (Xm'| ^comp^m'^comp |x m <) • (9) 

where Ck is the number of molecules which initially have been in the state k 
and m' - the index of an arbitrary chosen molecule in state k. The energy 
levels population numbers Ck are defined by suitable distribution function, for 
example, by the Boltzmann distribution, and satisfy the condition 

E k C k = M. (10) 
The expression (^) can be written in the form 

< J x >^ MSp(U+ mpPm ,U comp I^) (11) 
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Pm ' = (l/M)J2tiCk\x k m >) (X k m ,\ ■ (12) 

From formal point of view the operator ([l2]) can be interpreted as a statistical 
operator of the molecule m! ', which initially was in the mixed state of the kind 
of (g) with Pfe = Cfe/M. This fact tempts to analyze the states separability 
by considering the individual molecule ml with its "averaged" density matrix 
( p"2| ) and multiplying the result by the total molecules number M. However 
such an operation replaces quantum stochasticity with the classical one. As 
it was pointed out in the introduction, the molecule spin levels populations 
Cfe and therefore the probabilities with different k are approximately equal 
. Then formally density matrix (|l2|) belongs to the proximity of maximally 
mixed state. The misuse of the statistical ( "averaged" over sample) density 
matrix (|l^) for description of real quantum states of an individual molecule 
leads to the conclusion that real individual wave functions remain separable 
during computation. The important point is that the statistics, used in (|l2|), 
describes the number of molecules on the different energy levels, whereas the 
quantum mechanical state of an individual molecule remain pure state and can 
be transformed to entangled state following the specific computation algorithm. 

The work was supported by the NIOKR RT Foundation under the grant No. 
14-29. 
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